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Abstract 

We review novel methods to investigate, control and manipulate 
neutral atoms in optical lattices. These setups allow unprecedented 
quantum control over large numbers of atoms and thus are very promis- 
ing for applications in quantum information processing. After intro- 
ducing optical lattices we discuss the superfluid (SF) and Mott insulat- 
ing (MI) states of neutral atoms trapped in such lattices and investi- 
gate the SF-MI transition as recently observed experimentally. In the 
second part of the paper we give an overview of proposals for quan- 
tum information processing and show different ways to entangle the 
trapped atoms, in particular the usage of cold collisions and Rydberg 
atoms. Finally, we also briefly discuss the implementation of quantum 
simulators, entanglement enhanced atom interferometers, and ideas for 
robust quantum memory in optical lattices. 

1 Introduction 

The advent of experimental creation of Bose-Einstein condensates (BEC) 
in dilute weakly interacting gases of bosonic Alkali atoms in 1995 has brought 
new and exciting prospects to the fields of atomic physics and quantum state 
engineering. All the atoms in a BEC show identical quantum properties and 
in many respects a BEC therefore behaves like one 'big' quantum particle. 
It also can be imaged without being destroyed, its interactions with laser 
light are much stronger than that of a single particle and a BEC is less 
fragile than one might expect. Experimentally observed life times are of the 
order of several seconds; BEC's can be moved, shaken and rotated without 
immediately destroying their striking quantum features. Therefore they 
turned out to be an ideal testing bed for fundamental quantum physics, 
basic techniques of quantum state engineering and for investigating atomic 
properties [2] . The successful achievement of BEC and the early studies of 
their properties led to to the Nobel Prize in 2001 for the pioneers in this 
field E.A. Cornell, W. Ketterle, and C.E. Wieman 0. 
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Figure 1: Laser setup and resulting optical lattice configuration in 3D. 



The properties of a BEC described above make it an ideal tool for cold 
atom physics and quantum state engineering by quantum optical methods. 
One of the most promising experiments for future applications is the loading 
of a BEC into an optical lattice [1J03IIU- Optical lattices are periodic con- 
servative trapping potentials which are created by interference of travelling 
laser beams yielding standing laser waves in each direction (see Fig. The 
laser light induces an AC-Stark shifts in atoms and thus acts as a conserva- 
tive periodic potential. The usage of a BEC for loading has the advantage 
that its large atom density allows a filling of few particles per site n > 1 
while laser cooled atoms loaded into an optical lattice typically only allow a 
filling factor smaller than one. Furthermore the atoms loaded from a BEC 
are ultracold at temperatures very close to zero so that they practically be- 
have as if their temperature was T = 0, in particular all of them occupy the 
lowest Bloch band. The basic setup of optical lattices and the loading of a 
BEC will be described in detail in Sec. EJ 

An important novel feature which comes about because of the large filling 
factor n > 1 is that interatomic interactions - due to s-wave collisions - of 
two or more atoms occupying the same lattice site become important. The 
interatomic interaction potential, which usually leads to incoherent collisions 
in a thermal cloud of atoms and is responsible for the mean field in a BEC, 
causes a coherent energy shift U of two particles occupying a single lattice 
site. For the most frequently used alkali atoms 87 Rb and 23 Na the interaction 
is repulsive U > and thus the interaction energy competes with the kinetic 
energy J gained by a particle when hopping from one lattice site to the next. 
Both of these parameters depend on the laser intensity. When increasing the 
laser intensity the barriers between neighbouring sites increases in height and 
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J goes down. At the same time two particles sitting in one site are pushed 
together closer resulting in an increased interaction energy. Therefore the 
competition between kinetic and interaction energy can be controlled by the 
laser intensity which is an external easily changeable control parameter. In 
Sec. EH we will describe the physical consequences of different ratios J/U and 
investigate the properties of the resulting ground states. 

It turns out that the nature of the resulting quantum ground states 
allows for a number of applications ranging from the controlled creation 
of molecules [7], the study of strongly correlated quantum states jH], en- 
hanced atom interferometry experiments and precision measurements |27j . 
controlled entanglement creation, to quantum computing |Sl I§1 ITU1 1121 I13| . 
In addition to having well defined initial motional states a number of differ- 
ent internal (hyperfine) atomic states can be trapped which results in addi- 
tional well controlled degrees of freedom. Optical lattices are very versatile, 
easily altered by changing the laser setup and the control over the quan- 
tum state of the atoms by quantum optical techniques is unprecedented. 
Furthermore, decoherence times in optical lattices are large compared to 
typical experimental times which allows the study of Hamiltonian dynam- 
ics of strongly correlated quantum states with time dependent parameters 
and enables quantum information processing. All of these features make 
optical lattices a unique system to study. We will discuss some of the possi- 
ble applications in Sec. |I] focussing on applications in quantum information 
processing. We will show how to implement quantum memory and quan- 
tum gates in optical lattices and review ideas for using optical lattices as a 
quantum simulator for complex strongly correlated systems. 

The high degree of precision that can be obtained in manipulating quan- 
tum optical systems and the large degree of isolation from the environ- 
ment has led to a number of different quantum optical setups which might 
be suitable for quantum computing. Several experiments in the groups of 
D. Wineland 15 a and R. Blatt ^7J on ion traps have convincingly demon- 
strated two qubit gates. The next major step in these efforts is to scale them 
to a larger number of ions which will probably be on the order of ~ 30 ^0] • 
Also, photons in an optical cavity might be used to implement a quantum 
computer j!4j . The biggest advantage of cavity systems is that they might 
allow to couple resting qubits (atoms) to flying qubits (photons) and thus 
enable the implementation of a quantum network. 

There are also several non-quantum optical systems that have been iden- 
tified as candidates for implementing a quantum computer. At the mo- 
ment the most advanced system is NMR (Nuclear Magnetic Resonance) 
|18j . However, although NMR systems can be used to perform demonstra- 
tions of quantum computing algorithms on very small systems it is as yet 
not known how to scale an NMR quantum computer. Other promising can- 
didates are superconducting Josephson junction arrays ^H], quantum dots 
|2U| and somewhat more exotic systems like e.g. ammonia molecules con- 
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Figure 2: Interaction of a laser with a single atom. The red (blue) arrows 
indicate red (blue) detuned lasers. 



fined inside fullerenes [U- Which of these setups will eventually lead to 
a successful implementation remains to be seen. However, as we will show 
for optical lattice setups in the remainder of this paper a lot of interesting 
quantum physics is contained in the above systems. Therefore, even if some 
of them will not succeed in performing scalable quantum computations, it 
is certainly worth investigating them to learn more about their astonishing 
quantum features. 

2 Loading a BEC into an optical lattice 

In this section we first briefly describe the creation of an optical lattice; 
specifically we show how the interaction between atoms and a laser field can 
lead to a conservative trapping potential. Then we investigate how such a 
lattice can be used to trap and manipulate neutral atoms from a BEC. We 
will also mention some of the complications and imperfections that arise in 
this context but not deal with them in detail. 

2.1 Optical lattices 

Let us for simplicity consider a two-level atom with an internal (hyper- 
fine) ground state \g) and an excited state |e) in one spatial dimension 
only coupled to a monochromatic classical laser field with a detuning A as 
schematically shown in Fig. |2j The laser is assumed not to cause significant 
population in the excited state |e) and therefore we may neglect sponta- 
neous emission from it. In a frame rotating with the laser frequency the 
Hamiltonian (with h = 1) for the atom is given by 

tfatom = + A |e> < C | - ^ (\e) (g\ + \g) (e\) . (1) 

Here x and p are the coordinate and momentum operators respectively, and 
the mass of the atom is to. The laser drives the transition between the 
two atomic states with a Rabi frequency f2 which is proportional to the 
laser field and the dipole transition matrix element of the atom. We now 
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write the atomic wave function as |^(t)) = \^f e (t)) |e) + |\P g (t)) \g) where 
|^e(g)(i)) = ( e (9)\ ^(t)} and derive the equations of motion for the wave 
functions |* e (t)) and \^f g (t)). We find 

: '£r + A )\*e(t))--y-\%(t)). (3) 



\2r7i 

If the laser is far detuned, i.e. |A| 3> \Q\ it is permissible to adiabatically 
eliminate the excited internal state of the atom \e). The adiabatic elimina- 
tion consists in assuming that the wave function \ ^ e (t)) adiabatically follows 
any change in the ground state wave function \^f g (t)). We can thus set the 
left hand side of Eq. © equal to zero. Furthermore we assume the kinetic 
energy of an excited atom (and thus also the recoil energy E R = k 2 /2m 
where k = 2tt/X with A the laser wave length) to be much smaller than the 
detuning En <C |A| and neglect it in Eq. (J3J) compared to A. The equation 
Eq. (J2J) is now easily solvable and substituting its solution into Eq. (J2J) we 
find that the atom now moves according to the Schrodinger equation 



dt \ 2m 

where we have defined the optical potential as 

We will always assume the atom to be interacting with a standing light wave 
which arises from the interference of two counter propagating laser beams. 
This means that the spatial dependence of the Rabi frequency is given by 

Q(x) = r2osin(/cx), (6) 

where f^o is a constant depending on the laser intensity and the properties 
of the atom. We therefore find that the motion of the atom is effectively 
described by a periodic conservative potential with period a = A/2 and 
a depth controlled by the laser intensity. For a sufficiently deep lattice 
the atoms can be trapped near the potential minima which are also called 
lattice sites. In the vicinity of each lattice site the trap potential is well 
approximated by a harmonic oscillator potential with frequency 

LOT = , • (7) 

V2A^ 
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2.2 Imperfections 

Above we have totally neglected spontaneous emission from the excited 
atomic level. We give an estimate of this effect by assuming that the atom 
is localized close to the minima of the optical lattice potential. For blue 
detuned optical lattices (A < 0) the lattice sites are at positions of zero 
laser intensity while for a red detuned optical lattice (A > 0) they are at 
locations of maximum laser intensity. Therefore in a blue lattice the effects 
of spontaneous emission on trapped particles are less strong than for red 
lattices. In a blue lattice we find the effective spontaneous emission rate j^s 
of a trapped atom 

Teff ~ ~^T, (8) 

where 7 is the spontaneous emission rate from the excited atomic level. 
Therefore spontaneous emission will be negligible for large detunings A. 
Indeed, it has been shown in an experiment by Friebel et al. [221 that the 
spontaneous emission rate in far detuned optical lattices can be of the order 
of several minutes. 

Since we are interested in having a few atoms per lattice we have to 
consider the effects of light induced inelastic collisions between two atoms. 
In Fig. [3] we have schematically shown the interaction potentials for the 
different internal states (S + S) and (P + S) of two atoms approaching each 
other at a distance r. As a consequence of this interaction the detuning 
A of the laser field changes as the particles come close to each other and 
becomes A = at the so called Condon point Rq- The probability for the 
two particles to be excited to the states P + S increases around Rc- The 
atoms can then either form a quasistable molecule (for red detuning) which 
decays by spontaneously emitting a photon or they can be promoted to an 
unbound S + P state with a large kinetic energy equal to the detuning of 
the laser at r = 00. Both of these effects should be kept as small as possible 
to minimize loss from the optical lattice. There are several ways to do so: 

• Avoid a Condon point by detuning far to the red. However, this will 
lead to a larger spontaneous emission rate compared to the case of a 
blue lattice. 

• Avoid a red detuning that puts the Condon point close to a quasi 
bound molecular state (shown in Fig. |2J) as this enhances the proba- 
bility of creating such molecules. 

• Choose the detuning such that the two potential curves have the 
largest possible angle with each other at the Condon point. This makes 
the probability of finding the two particles in the region where A ~ 
small. 

The above conditions for neglecting spontaneous emission and making laser 
induced collisions unimportant can be fulfilled quite well experimentally and 
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Figure 3: Schematic interaction potential curves for two atoms at distance r. 
The spontaneous emission from a bound molecular state 7 is also indicated. 



then our above simple model for an optical lattice as a conservative potential 
is valid. 

2.3 Experimental setups 

The situation where one pair of counter propagating laser beams interferes 
and yields a trap potential in ID while the atoms are only weakly trapped 
by a magnetic field (or not trapped at all) in the other two dimensions is 
usually called a ID optical lattice. A BEC was first loaded into a ID lattice 
by B.P. Anderson and M.A. Kasevich [3]. In this experiment a magnetically 
trapped BEC is loaded into a weak vertical optical lattice while at the same 
time turning off the magnetic trapping almost completely. The experiment 
is designed such that the interplay between the gravitational force and the 
optical lattice leads to a peculiar interference pattern of atoms falling out 
of the lattice which allows to study the coherence properties of the original 
BEC. 
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Figure 4: Experimental setup used in Munich to load a BEC into a 3D 
optical lattice [Hj. The BEC is produced in a magnetic trap located at the 
left hand side of the figure and is then transported to the right where optical 
access from all directions is possible as shown by the red arrows. 



A ID optical lattice setup is theoretically easily extended to the 3D 
situation as shown in Fig. ^ Three pairs of counter propagating laser beams 
interfere to produce an optical potential which traps the atoms in each 
dimension. Experimentally, however, the achievement of this situation is 
less straightforward. The experimental setup has to allow for optical access 
to the BEC from three different orthogonal directions and this is usually 
incompatible with producing the magnetic trap required to create the BEC. 
For this reason the group headed by T. Hansch in Munich realized the 
concept of spatially separating the production of a BEC and performing 
the actual optical lattice experiment. Their setup is schematically shown in 
Fig. 0] A number of coils are used to create an inhomogeneous magnetic 
field that guides the BEC out of its original position and transports it into 
the 'science chamber' where optical access from three directions is possible 
and allows the realization of a 3D optical lattice. Next we will discuss the 
novel features of a BEC loaded into such a 3D optical lattice. In particular 
we will investigate the effects of interactions between the particles which are 
responsible for striking new features in comparison to previous experiments 
on optical lattices. 
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3 Superfluids and Mott insulators 



The above setup with 3D lattice geometry was predicted to yield a qual- 
itatively new behaviour of a dense cloud of ultracold atoms in 1998 |23j 
and was experimentally realized in 2001 by M. Greiner et al. We can 
gain some insight into the properties of the atoms trapped in the lattice by 
looking at two limiting cases; the superfluid limit where the kinetic energy 
dominates the repulsive interaction and the opposite limit which is called 
the Mott insulating limit. 

3.1 Superfluid limit 

When the optical lattice is shallow the kinetic energy dominates over the 
repulsion of two particles sitting in the same lattice site. In this situation 
nearly all the atoms occupy the same single particle state. For very small 
interaction energies this single particle state is very close to the Bloch wave 
function with quasi momentum q = of the lowest Bloch band, i.e. the single 
particle ground state of the periodic potential. As this wave function is de- 
localized over the system there is long range off-diagonal coherence over the 
whole lattice. States like this are very well described by the Gross-Pitaevskii 
equation 2_ which takes the interaction into account via a simple mean field 
approximation. Neglecting interactions altogether we can understand the 
form of the underlying Hamiltonian intuitively. We define a bosonic de- 
struction operator ai destroying an atom in the lattice site labelled I. Its 
adjoint operator a\ correspondingly creates a particle in this lattice site and 
n\ = a\cn counts the number of particles in site I. By hopping (tunneling) 
from one lattice site to an adjacent site an atom loses a kinetic energy J de- 
pending on the height of the barrier between the sites and the Hamiltonian 
therefore reads 



where the sum runs over all pairs of nearest neighbours in the lattice. An 
atom in the lattice will therefore minimize its energy when it is in a super- 
position state of many lattice sites. 

A more detailed investigation shows that the state 



of N particles with |vac) the vacuum minimizes the energy. This state clearly 
corresponds to each particle occupying the lowest Bloch state and spreads 
the wave function of each particle over the whole lattice as indicated in 
Fig. EH. The state therefore shows large off diagonal long range coherence 
as can be seen by correlating the destruction of an atom in lattice site I with 
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Figure 5: Schematic representation of a) a SF state and b) a MI state with 
commensurate filling n = 1 in an optical lattice with M sites. 



the creation of a particle in site m. The corresponding correlation function 
is found by the expectation value p; >m = (a]a m ). In the SF regime this 
correlation function is independent of the distance between sites I and m and 
its modulus is given by \pi t m\ = N/M where M is the total number of lattice 
sites. We also note that the particle number fluctuations (An;) 2 = N/M in 
this state are large. 

The lowest lying excitations in the SF limit are determined by the Bloch 
band and consist of waves with quasi momentum q. The energy of the 
first excited state scales like t\ oc J/M 2 which tends to for large systems 
M — ► oo. This vanishing gap between the ground state and low lying excited 
states makes the SF state rather fragile and seems to render it very difficult 
to adiabatically load a BEC into an optical lattice. Fortunately, however, 
the turning on of a lattice does not induce quasi-momenta and therefore does 
not excite low lying excitations. Therefore loading a lattice from a BEC is 
possible at reasonable speed without significant heating of the atomic cloud. 

3.2 Mott insulating limit 

By increasing the lattice depth the kinetic energy goes down because the 
barrier height increases. At the same time the size of each lattice site de- 
creases which causes the atoms to repel each other more strongly. Again we 
can understand this situation intuitively. We neglect any kinetic energy and 
only consider the interaction energy of two particles in the same lattice site. 
As usual the two particle interaction energy increases with the square of the 
particle density which itself is proportional to the square of the number of 
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particles per well. The Hamiltonian therefore is given by 

h mi = ^ S>,(n,-1). (11) 

i 

Here the parameter U is the interaction energy of two particles in the same 
lattice site. Note the —1 term which is usually neglected in most mean 
field theories. It ensures that a single particle does not gain any interaction 
energy. 

In this situation particles in the optical lattice minimize their energy by 
staying away from each other instead of spreading out over the whole lattice. 
Let us consider the case of commensurate filling M = N where the smallest 
interaction energy is obtained by putting each particle into its individual 
lattice site. A state like this is written as 

i^mi) «n a ii vac )' ( i2 ) 
i 

is called a Mott insulator, and schematically shown in Fig. It does not 
have any interaction energy (\&mi| ^int I^mi) = while the above SF state 
would have an interaction energy of (^sfI -f^int |^sf) = UN(N — 1)/2M 2 . 
Furthermore, the MI state has no long range coherence pi m = for I ^ m 
and also the particle number fluctuations are zero Anf = 0. 

The MI state is much more stable than the SF state. The lowest lying 
excitations are obtained by taking one particle out of its site and putting it 
on top of an atom sitting in another lattice site. The energy required for 
this is U independent of the size M , i.e. the system is characterized by a 
finite gap U large compared to the small excitation energies in the SF limit. 

Experimentally the situation of commensurate filling can be achieved 
by suitably superimposing the optical lattice with a much weaker magnetic 
trapping potential. The atoms fill those lattice sites which are in the center 
of the wide trap and avoid the outer ones so that we obtain commensurate 
filling in the center of the lattice. The magnetic trap acts as a 'soft' boundary 
which adjusts to the number of particles. Note, however, that the nature of 
the MI will be altered close to this boundary. 

3.3 Intermediate regime 

If the two parameters J and U are of similar magnitude the situation is 
more complicated. Kinetic energy and interactions compete with each other 
and only a detailed investigation shows whether the system is a SF or a MI. 
For commensurate filling of the lattice with a filling factor of n = 1 one 
easily finds that the interaction wins for U > U c while the kinetic energy 
dominates for U < U c , where the critical point is U c ~ 5.8zJ. Here z = 2d 
is the number of nearest neighbours in a d-dimensional lattice. For larger 
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Figure 6: Phase diagram of the 3D Bose-Hubbard model showing the MI 
regions with filling factor n and the SF region for different ratios J/U and 
chemical potentials fi/U. Note that the simple mean field approximation 
used here does not yield U c = 5.8zJ. 



commensurate fillings of the lattice the value of U c increases because the 
interaction energy goes up quadratically with the filling while the gained 
kinetic energy only increases linearly with the filling factor. 

For incommensurate fillings the system is best described in the grand 
canonical ensemble where the mean particle number is fixed by a particle 
reservoir of chemical potential [i. By using simple mean field theory the 
SF and MI regions can be found approximately; the numerical result is 
shown in Fig. When varying J over the critical point the ground state 
abruptly changes its nature and the transition between the SF and the MI 
regime at temperature T = takes place in a so called quantum phase 
transition if the lattice is sufficiently large. We can understand the basic 
physics of a quantum phase transition by considering a Hamiltonian of the 
form H = H\ + gUi where the two parts of H commute, i.e. [Hi, H2] = 
and g > is a parameter. Since Hi and H2 commute the eigenstates of H 
are independent of g while the eigenvalues change with g. For small values 
of g < g c the ground state of H will be determined by the ground state of 
Hi while for large values of g > g c it is given by the ground state of H2 ■ If 
these two ground states have totally different properties the system alters its 
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Figure 7: Interference patterns for decreasing J/J7 (a-h). For U < U c the 
long range coherence leads to an interference pattern which clearly indicates 
the band structure (a-d). The interference fringes vanish when the system 
enters the MI regime U > U c (e-h). 



behaviour suddenly when the parameter g is changed over the critical point 
g c . A quantum phase transition also occurs when the commutator of the 
two parts of the Hamiltonian is not exactly zero but does not increase with 
the size of the system and therefore can be neglected in the case of large 
systems. This is exactly the situation we encounter in the so called Bose- 
Hubbard model with Hamiltonian -Hbhm = #SF + #mi we are considering 
here for a changing ratio of J/U . 

An optical lattice has the unique feature of permitting studies of the 
dynamics of the quantum phase transition in Hbkm, i-e., the parameters 
U and J can be changed on a time scale much smaller than the coherence 
time of the system but still sufficiently slow to avoid significant excitations. 
The experiments carried out by M. Greiner et al. |S] to study this dynamics 
started out from a Bose-condensate in a magnetic trap. Then the optical 
lattice was (adiabatically) turned on. Its depth determines the ratio J/U and 
thus the kind of ground state in the optical lattice. To measure the properties 
of the atomic state the lattice potential is then turned off and the particles 
are allowed to expand freely before their density distribution is measured. 
In the case of a SF state which has long range off diagonal coherence the 
particles interfere and the resulting interference pattern clearly indicates the 
band structure of the lattice. When the system is in a MI state there is no 
long range coherence and therefore the interference pattern disappears. This 
behaviour is similar to interference in optics experiments where interference 
patterns also vanish when there is lack of coherence. The corresponding 
interference patterns as measured in 0| are shown in Fig. for different 
depths of the optical lattice. 

Next we discuss some of the applications of MI states in optical lattices 
focussing on applications in quantum computing. Afterwards we look at 
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atomic interferometry enhanced by using multi particle entangled states. 

4 Quantum information processing in optical lat- 
tices 

The experiments on the SF to MI transition are a very important step to- 
wards realizing theoretical ideas for controlled entanglement creation in op- 
tical lattices by interactions between neutral atoms [51 11UI 151 112j. Together 
with single particle manipulations such entanglement operations constitute 
the basic building blocks for a quantum computer which fulfils all the 
basic requirements as set out by D. DiVincenzo [2^]. First important exper- 
imental steps towards realizing multi particle entangled states have recently 
been achieved by Mandel et al. 26 j. In this section we will first discuss 
the requirements for implementing a quantum computer and then show how 
they can be satisfied in an optical lattice. 

4.1 What is quantum information processing? 

A computation can be considered as a physical process that transforms an 
input into an output. A classical computation is that in which the physi- 
cal process is based on classical laws while the physical process in quantum 
computation is based on quantum laws and in particular on the superposi- 
tion principle. In quantum information processing the inputs and outputs 
are encoded into the states of the system and the computation evolves the 
input state by a designed well controlled unitary time evolution. In general, 
if our algorithm consists of evaluating a given function /, we can design an 
interaction Hamiltonian such that the resulting unitary evolution transforms 
the input states according to \ j) (g) |0) — > \j) |/(j)) (j = 1, N). 

Using a quantum computer we can do more than with a classical com- 
puter. When preparing the input state in a superposition of all L = 2 N 
basis states for N qubits it will be transformed according to 



in a single run. All the values of the function / are present in this superposi- 
tion. However, we do not have access to all of this information since after a 
measurement we will only obtain one result with a certain probability. The 
property of using quantum superpositions to run the computer only once 
was termed quantum parallelism by Feynman. 

4.1.1 Requirements 

A quantum computer consists of a quantum register (in our case an array 
of atoms in the optical lattice) that can be manipulated and measured in 
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a controlled way. If we want to build a quantum computer, we need the 
following elements 25_: 

(i) A set of qubits which are two-level systems forming the quantum 
register. In optical lattices this set of qubits is formed by two internal states 
of the atoms trapped in the lattice. We denote by {|0) ■ , 1 1) •} two orthogonal 
states of the j-th qubit (also called the computational basis), so that a state 
\^>) of all the qubits in the quantum register can be written as 

i 

I*) = c hh-iN\h)i®\h)2® \Jn) n 

jl J2vJJV=0 
1 

= Yl C hh-3N \h,h,---jN) • (14) 

jl J2vJJV=0 

These qubits can be in superposition and entangled states, which gives the 
extraordinary power to the quantum computer; they also have to be well 
isolated from the environment to prevent decoherence. 

(ii) A universal set of quantum gates which allows the controlled ma- 
nipulation of the qubits according to any unitary operation U so that 
1^) — * U\^). Fortunately, this task is enormously simplified given the fact 
that any U can be decomposed as a product of gates belonging to a small 
set, a so-called universal set of gates, i.e. if we are able to perform the gates 
of this set we will be able to perform any unitary operation on the register 
by applying a sequence of them. There are many sets of universal gates 
and we will concentrate on the set containing one two qubit phase gate and 
single qubit operations as this is the one most naturally arising in optical 
lattice setups. Note that two-qubit gates require interactions between the 
qubits, and therefore are the more difficult ones in practice. 

(iii) Detection of the output state. One should be able to measure each 
qubit in its computational basis. This process requires the interaction with 
a measurement apparatus in an irreversible way and is thus not governed 
by a unitary time evolution. 

(iv) Initial state preparation. We must be able to erase quantum registers 
to prepare the initial quantum state, for example the state |0, 0, . . . ,0). This 
is actually not an extra requirement, it is enough to detect the qubits and 
to apply a single-qubit gate flipping the qubit |1) — > |0) and |0) — > |1) if 
necessary. 

(v) Scalability of the system. The difficulty of performing gates, mea- 
surements, etc., should not grow (exponentially) with the number of qubits. 
Otherwise, the gain in the quantum algorithms would be lost. 

(vi) In addition to the above we also require networking ability. The 
static qubits used to perform quantum computations should be trans- 
formable into flying qubits which can easily be transmitted between different 
locations. 
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In the following we will discuss how to implement an initialized quantum 
register and a universal set of quantum gates in an optical lattice setup. 
Due to its periodic structure the optical lattice provides an intrinsically 
scalable setup. The problem of detecting the output state is closely related 
to addressing single qubits as required for performing single qubit operations. 
A discussion of how to achieve networking ability lies beyond the scope of 
this paper. 

4.2 Quantum memory and single qubit gates in optical lat- 



in optical lattices two states of the ground state hyperfine manifold of the 
trapped atoms \a) = |0) and \b) = |1) are ideal candidates to implement a 
qubit. However, in order to be able to manipulate and control atomic qubits 
it is necessary to know the position of each of the atoms precisely. While 
random filling of optical lattices from laser cooled atoms and a superfluid 
filling of a lattice (with large particle number fluctuations as discussed above) 
do not provide sufficient knowledge on the position of the atoms a MI state is 
ideally suited for this purpose. All the lattice sites are filled, i.e. each lattice 
site contains a qubit and the fluctuations in the occupation numbers are 
very small. For an appropriate choice of internal atomic states the optical 
lattice allows to trap both of them and thus can hold a quantum register for 
storing quantum information. 

In principle it is straightforward to induce single qubit gates by using 
Raman transitions between the two internal states \a) and \b). Raman tran- 
sitions with Rabi frequency and detuning 5 are described by the Hamil- 
tonian 



which induce rotations of the qubit state on the Bloch sphere. The axis 
and angle of this rotation depend on the choice of laser parameters and can 
be chosen freely. The major problem in inducing single qubit operations 
is the addressing of a single atom as it is difficult to focus a laser to spots 
of order of an optical wave length which is the typical separation between 
atoms in the lattice. Possible solutions to this difficulty are using schemes 
for pattern loading (2211201) where only specific lattice sites are filled with 
atoms (cf. Fig. [SJ) or using additional marker atoms which specify the atom 
the laser is supposed to interact with. 

MI atoms in an optical lattice have already experimentally been used 
as qubits and it has been shown that they support multi particle entangled 
states [2H1 • The robustness of these qubits is, however, limited by stray 
magnetic fields and spin echo techniques need to be used to perform the 
experiments successfully. A different approach to obtaining robust quantum 
memory uses a more sophisticated encoding of the qubits .271. A ID chain 
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Figure 8: Selective filling of optical lattice sites in a superlattice structure. 
The red boxes indicate lattice sites filled with exactly one particle. The blue 
boxes show lattices sites that are empty. 



with an even number of atoms encodes a single qubit in the states |0) = 
\abab---ab) and |1) = \baba---ba). These states both contain the same 
number of atoms in internal states \a) and \b) and therefore interact with 
magnetic fields in exactly the same way avoiding dephasing of the quantum 
information stored in the chain. This and the fact that all the atoms have 
to flip their internal state to get from one logical state to the other make 
these qubits very robust against the most dominant experimental sources 
of decoherence. However, at the same time this robustness makes it more 
difficult to manipulate them. A laser pulse no longer corresponds to a simple 
rotation on the Bloch sphere spanned by the logical states |0) and |1) which 
makes the realization of a single qubit gate difficult. We will describe how 
these problems can be circumvented in Sec. 14.51 and further details can be 
found in |27j . 
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4.3 Two qubit operations 

Implementing a two qubit gate is more challenging than the single qubit 
gates. The different schemes for two qubit gates can be classified in two 
categories. The first version relies on the concept of a quantum data bus; 
the qubits are coupled to a collective auxiliary quantum mode, like e.g. a 
phonon mode in an ion trap, and entanglement is achieved by swapping the 
qubits to excitations of the collective mode. The second concept which is the 
basis for two qubit gates between atoms in optical lattices deploys control- 
lable internal-state dependent two-body interactions. Examples for different 
interactions are coherent cold collisions of atoms, optical dipole-dipole inter- 
actions [HI El and the 'fast' two-qubit gate based on large permanent dipole 
interactions between laser excited Rydberg atoms in static electric fields jTUj . 
Besides these dynamical schemes for entanglement creation it is also possible 
to generate entanglement by purely geometrical means [2B]- We will now 
discuss the different ways to achieve two qubit operations in optical lattices 

4.3.1 Entanglement via coherent ground state collisions 

The interaction terms introduced in Eq. (|11|) which describe s-wave col- 
lisions between ultracold atoms in one lattice site are analogous to Kerr 
nonlinearities between photons in quantum optics. For atoms stored in op- 
tical lattices these nonlinear atom-atom interactions can be large jSj, even for 
interactions between individual pairs of atoms, thus providing the necessary 
ingredients to implement two-qubit gates. 

We consider a situation where two atoms in a superposition of internal 
states | a) and \b) are trapped in the ground states of two optical lattice 
sites (see Fig. Initially, at time t = — r these wells are centered at 

positions sufficiently far apart so that the particles do not interact. The 
optical lattice potential is then moved state selectively and for simplicity 
we assume that only the potential for a particle in internal state \a) moves 
to the right and drags along an atom in state \a) while a particle in state 
| b) remains at rest. Thus the wave function of each atom splits up in space 
according to the internal superposition of states \a) and \b). When the wave 
function of the left atom in state \a) reaches the second atom in state \b) 
as shown in Fig. they will interact with each other. However, any other 
combination of internal states will not interact and therefore this collision 
is conditional on the internal state. A specific laser configuration achieving 
this state dependent atom transport has been analyzed in Ref. [H] for Alkali 
atoms, based on tuning the laser between the fine structure excited p-states. 
The trapping potentials can be moved by changing the laser parameters. 
Such trapping potentials could also be realized with magnetic and electric 
microtraps |29] , 

We therefore only need to consider the situation where atom 1 is in state 
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atom 1 



atom 2 



a) 




x a (t) 



b) 




ve-v 



z 6 (t)=const 



Figure 9: We collide one atom in internal state \a) (filled circle, potential in- 
dicated by solid curve) with a second atom in state \b) (open circle, potential 
indicated by dashed curve)). In the collision the wave function accumulates 
a phase according to Eq. (fTT)) . a) Configurations at times t = ±r and b) at 
time t. 



| a) and particle 2 is in state \b) to analyze the interactions between the two 
atoms. The positions of the potentials are moved along trajectories x a (t) 
and x b (t) = const so that the wave packets of the atoms overlap for a certain 
time, until they are finally restored to the initial position at the final time 
t = t. This situation is described by the Hamiltonian 



and u is the atom-atom interaction term (which lead to the interaction 
term in Eq. Qll|) )- Ideally, we want to implement the transformation from 
before to after the collision, 



r (x a - x a (-r))4(x b - x\-r)) -> e^ a o(x a ~ x a (r))4(x b - x\r)), (17) 



where each atom remains in the ground state ipQ of its trapping potential 
and preserves its internal state. The phase cj) = <ft a + <p b + ab will contain a 
contribution </> ab from the interaction (collision) and (trivial) single particle 
kinematic phases (j) a and <fi b . The transformation Eq. (|17|) can be realized 
in the adiabatic limit, whereby we move the potentials slowly on the scale 
given by the trap frequency, so that the atoms remain in their motional 




/3=a,b L 



(16) 




are position and momen- 
describe the displaced trap potentials 



19 



ground state. In this case the collisional phase shift is given by (j) = 
dtAE(t)/h, where AE(t) is the energy shift induced by the atom-atom 
interactions 

AE(t) = J dx\r a (x - x a (t)) \ 2 \4 (x - x b (t)) | 2 , (18) 

with a s the s-wave scattering length. In addition we assume that |A.E(i)| <C 
Hujt so that no sloshing motion is excited. 

The interaction phase thus only applies when the atoms are in internal 
state | a, b) but not otherwise. Carrying out the above state selective collision 
with a phase cff" h = it we obtain (up to trivial phases) the mapping (as above 
we identify \a) = |0) and \b) = |1)) 

|0,0) -> |0,0) 

|0,1) - -|0,1> 

1 1, 0> -» 1 1, 0> 

|1,1) -> |1,1) (19) 

which realizes a two-qubit phase gate that is universal in combination with 
single-qubit rotations. 



4.3.2 State-selective interaction potential 

An alternative possibility, for a nontrivial logical phase to be obtained, is to 
rely on a state-independent trapping potential, while defining a procedure 
where different logical states couple to each other with different energies. 
An example is given by the interaction between state-selectively switched 
electrical dipoles [TU] . 

In each qubit, the hyperfine ground states \a) = |0) is coupled by a laser 
to a given Stark eigenstate \r) which does not correspond to a logical state. 
The internal dynamics is described by a model Hamiltonian 



iT/(i,xi,x 2 ) = ^2 

i=i,2 



5Mr)j(r\-^^(\a) j (r\ + h.c.) 



2 

+ «|r)i(r| <8) \r) 2 (r\, (20) 



with Qj(t, Xj) Rabi frequencies, and Sj(t) detunings of the exciting lasers. 
Here, u is the dipole-dipole interaction energy between the two particles. We 
have neglected any loss from the excited states \r)j. We discuss two possible 
realizations of two qubit gates with this dynamics. The most straightforward 
way to implement a two-qubit gate is to just switch on the dipole-dipole 
interaction by exciting each qubit to the auxiliary state |r), conditioned on 
the initial logical state. This can be obtained by two resonant (5% = 62 = 0) 
laser fields of the same intensity, corresponding to a Rabi frequency f2i = 
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0,2 ~S> u. After a time r = p/u, the gate phase <p is accumulated and the 
particles can be taken again to the initial internal state. However, besides tp 
being sensitive to the atomic distance via the energy shift u, during the gate 
operation (i.e. when the the state \rr) is occupied) there are large mechanical 
effects, due to the dipole-dipole force, which create unwanted entanglement 
between the internal and the external degrees of freedom. These problems 
can be overcome by assuming single-qubit addressability and by moving to 
the opposite regime of small Rabi frequencies £li(t) ^ ^(i) -C u. The gate 
operation is then performed in three steps, by applying: (i) a 7r-pulse to 
the first atom, (ii) a 27r-pulse (in terms of the unperturbed states) to the 
second atom, and, finally, (hi) a 7r-pulse to the first atom. The state 1 00) 
is not affected by the laser pulses. If the system is initially in one of the 
states 1 01) or 1 10) the pulse sequence (i)-(iii) will cause a sign change in 
the wave function. If the system is initially in the state |11) the first pulse 
will bring the system to the state i\rl), the second pulse will be detuned 
from the state |rr) by the interaction strength u, and thus accumulate a 
small phase <p ~ irQ2/2u *C 7r. The third pulse returns the system to the 
state e l ( n ~^\ll) , which realizes a phase gate with (p = tt — (p ~ tt (up to 
trivial single qubit phases). The time needed to perform the gate operation 
is of the order r ~ 2tt/VL\ + 2tt/£12- Loss from the excited states \r)j is 
small provided "fAt <C 1, i.e. IX,- 3> 7. A further improvement is possible 
by adopting chirped laser pulses with detunings £1,2 (*) = S(t) and adiabatic 
pulses ^1,2^) = Q(t), i-e., with a time variation slow on the time scale given 
by Q and 5 (but still larger than the trap oscillation frequency), so that the 
system adiabatically follows the dressed states of the Hamiltonian Hj. As 
found in ^U], in this adiabatic scheme the gate phase is 



tp(r) 



[ t0+T dt sgn(£) ^ ^l +gg! _ sgn (£) f\ S \- y/^T^ 

Jt 2 V 



(21) 



with 5 = 5 — 2 /(4i5 + 2u) the detuning including a Stark shift. For a spe- 
cific choice of pulse duration and shape Q(t) and 5(t) we achieve </?(r) = tt. 
To satisfy the adiabatic condition, the gate operation time r has to be 
approximately one order of magnitude longer than in the other scheme dis- 
cussed above. In the ideal limit Qj <C u, the dipole-dipole interaction energy 
shifts the doubly excited state |rr) away from resonance. In such a "dipole- 
blockade" regime, this state is therefore never populated during gate opera- 
tion. Hence, the mechanical effects due to atom-atom interaction are greatly 
suppressed. Furthermore, this version of the gate is only weakly sensitive to 
the exact distance between the atoms, since the distance-dependent part of 
the entanglement phase (p <C vr. For the same reason, possible excitations 
in the particles' motion do not alter significantly the gate phase, leading to 
a very weak temperature dependence of the fidelity. 

In principle the controlled implementation of the above single and two- 
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qubit quantum gates in the scalable setup of an optical lattice - avoiding 
any decoherence - allows any quantum computation to be performed, and 
in particular to run Shor's algorithm |M1| for factoring large numbers in 
polynomial time. We note that the implementation of this factoring algo- 
rithm is one of the major practical goals in quantum computing as most 
current cryptographic systems rely on the difficulty of this task. However, 
for quantum computers to exceed classical ones in tasks like factoring they 
would have to operate thousands of qubits and the technology for achieving 
this goal might only be reachable in years from now. Still, there are also 
important nontrivial short term goals which can be realized with nowadays 
technology for quantum information processing. We present one of them in 
the following section. 

4.4 Feynman's universal quantum simulator 

A universal quantum simulator (QS) is a controlled device which efficiently 
reproduces the dynamics of any other many-particle system that evolves 
according to short range interactions. Therefore a QS could be used to ef- 
ficiently simulate the dynamics of a generic many-body system, and in this 
way function as a fundamental tool for research in many body physics. Ex- 
amples include spin systems, where classical simulations become nontrivial 
for a few tens of spins, and also interacting electrons moving in a magnetic 
field in 2D, which show interesting effects like the Fractional Quantum Hall 
effect. The simulation of such nontrivial systems is well suited for atoms 
trapped in an optical lattice and might therefore be implementable with 
resources available in the lab in the near future. 

The realization of 2D motion in a 3D optical lattice is straightforward. 
Since the height of the barriers between sites adjacent in z-direction in- 
creases with the intensity of the corresponding pair of laser beams one can 
turn off any motion in z-direction by using high intensities for these lasers 
while keeping the barrier between sites adjacent in x, y direction low. Atoms 
will then only be able to tunnel along the x and y direction and will thus 
be confined to one layer of the optical lattice. The physical simulation of 
a magnetic field in z direction is less simple since time reversal symmetry 
needs to be broken. However, by accelerating the lattice and using addi- 
tional laser beams as shown in |34| the simulation of a magnetic field can be 
achieved and allows to study 2D atom gases in effective magnetic fields. As 
was shown by Hofstaedter [33] f° r strong magnetic fields the energy bands 
form a fractal structure depending on the parameter a = ABe/27rh, where 
A is the area of one of the elementary cells of the lattice, B is the strength of 
the magnetic field, and e is the charge of the particles. Plotting the eigenen- 
ergies as a function of a yields the so called Hofstaedter butterfly shown in 
Fig. 1101 While for electrons moving on a 2D surface it is very difficult to 
reach sufficiently large values of a to see this butterfly experimentally the 
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Figure 10: The Hofstaedter butterfly. Each dot in this fractal indicates an 
eigenenergy of the atoms moving in a magnetic field with parameter a as 
given in the text. The energy is normalized to eo which is given by the 
modulus of the hopping energy between two sites. 



parameter a can be varied easily between and 1 in an optical lattice and 
thus it should be possible to observe the energy band structures shown in 
Fig. EH 

The first idea for simulating (anti)ferromagnetic spin systems |32| in op- 
tical lattices came about after realizing that the interactions between atoms 
in optical lattices can be controlled [HI El - In a more general approach it was 
recognized by Jane et al. jllj that the nature of the controlled Hamiltonians 
available in optical lattices (and other quantum optical systems) is very well 
suited for simulating the time evolution of generic spin chains which are 
built up by qubits and evolve according to a Hamiltonian 

H N =j2H {i) +J2 Hm ( 22 ) 

that decomposes into one-qubit terms H^' and two-qubit terms H^ J ' . It 
can be shown that the time evolution resulting from Hn is decomposable 
into a series of short pieces K n which are generated by single- and two-qubit 
evolutions. The single qubit evolution is realized in an optical lattice by 
properly shining a laser on the trapped atoms and two-qubit evolutions are 
obtained by operating two-qubit quantum gates via atom-atom interactions. 
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A more detailed investigation shows that interesting spin-chain Hamilto- 
nians like the (anti) ferromagnetic Heisenberg Hamiltonian can be simulated 
without the need for local addressing of single atoms and could thus be 
realized with standard optical lattice setups. The ability to perform inde- 
pendent operations on each of the qubits would allow all possible bipartite 
Hamiltonians to be simulated and in particular the study of quantum phase 
transitions in spin chains. We will next briefly study an optical lattice setup 
with Rydberg interactions between adjacent atoms which shows a quantum 
phase transition and could be viewed as a quantum simulator [27]. We will 
also show how it can be used to implement an atom interferometer whose 
sensitivity is enhanced by the created entanglement. Furthermore, the ideas 
presented in the next section are the basis for implementing single qubit 
gates on robust qubits as introduced in Sec. 14.21 

4.5 Multi particle maximally entangled states in optical lat- 
tices 

Let us assume that the Rydberg interactions between neighbouring atoms 
are turned on all the time and that the atoms are arranged in a ID chain. 
Furthermore we allow each atom to either tunnel between two adjacent wells 
whose ground states are denoted by \a) and \ b) or to have an additional laser 
driving Raman transitions between two internal states (again denoted by \a) 
and 1 6)). This situation is schematically shown in Fig. 1111 

The dynamics of this setup can be understood as follows. Hopping be- 
tween the two modes |a) ■ and \b)j of the j-th atom is described by 

H h = B £ (V>, (b\ + h.c.) = B £ 4 j) (23) 

3 3 

with Ox the Pauli x-matrix for the two states (viewed as a spin) of the 
j-th atom and B > the energy associated with this process. The ground 
state of this Hamiltonian is obtained by putting each atom into an equal 
superposition of states written in spin notation as \[ x ) j = (If*) — || z ))/v / 2, 
where = \a) and || 2 ) = \b) i.e. we find 

I*/,) = \LLL •■■ ix) (24) 

Note that this part of the Hamiltonian is identical to a spin chain in a 
magnetic field B along the x-axis. 

The interaction between two nearest neighbours due to Rydberg inter- 
actions is given by 

H i = Wj2 (T z ) ^ +1) (25) 

3 

accounting for the energy difference of having two adjacent particles in the 
same vs. different internal states. In this case the ground state depends on 
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Figure 11: Creation of robust multi particle entangled states in ID beam 
splitter setups. Closed circles indicate an atom open circles an empty posi- 
tion. 



the sign of the interaction parameter W. For repulsive interactions W > 
the interaction will be minimized by arranging the particles as 



Both of these states are maximally entangled multi particle states, i.e. ex- 
tensions of the well known GHZ states for three particles. 

In the total Hamiltonian H = + Hi the interaction energy and the 
hopping energy compete with each other resulting in a quantum phase tran- 
sition. When the interaction energy is kept constant and the hopping term 
is switched off adiabatically, as e.g. shown in Fig.^^the state of the system 
will dynamically change from \^h) which is a product state to one of the 
two states \^i r ) or \^i a ) depending on the sign of the interaction. The exact 
values of the parameters a and (5 depend on the details of the dynamics 
and are discussed in [27j. These maximally entangled states can serve dif- 
ferent purposes depending on the sign of W. Let us discuss the possible 
applications of these two kinds of maximally entangled states 

4.6 Repulsive interactions 

In this case both parts of the superposition state have the same number 
of particles in each of the two internal states and thus external stray fields 
act identically on both, therefore not affecting the parameters a and (3. 
Because of this stability one can use them for storing quantum information 
in a robust way as described in Sec. 14.21 Single qubit gates can be performed 
by dynamically going back and forth through the quantum phase transition 
in the whole chain changing the parameters a and (3 in a controlled way. 




(26) 



while for attractive interactions the ground state is 




(27) 
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4.7 Attractive interactions 



For attractive interactions the terms in the superposition of \ipi a ) will re- 
spond to external fields very differently. Therefore the relative phase be- 
tween the parameters a and (3 will be very susceptible to these fields, in fact 
for N particles in the chain this phase will be N times larger than if there 
was just a single particle. The two parts of \ipi a ) can thus be used as two 
arms of an entanglement enhanced atomic interferometer. 



5 Summary 

This article reviewed recent advances in controlling and manipulating neu- 
tral atoms by quantum optical methods concentrating on the use of optical 
lattices. We investigated novel physical features arising from interactions 
in ultracold dense clouds of atoms as provided by BEC and showed how 
they enable experimentalists to gain control over the spatial degrees of free- 
dom of individual atoms. We also studied possible new applications - in 
particular in quantum information processing - which make use of this un- 
precedented quantum control. We showed that universal sets of quantum 
gates can be implemented and how they might lead to a universal scalable 
quantum computer in the future. Furthermore we also investigated possibil- 
ities for shorter term applications of optical lattices as quantum simulators 
for strongly correlated systems. 
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